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Abstract. We consider the behaviour of a critical system in the presence of a gradient 
perturbation of the couplings. In the direction of the gradient an interface region 
separates the ordered phase from the disordered one. We develop a scaling theory for 
the density profiles induced by the gradient perturbation which involves a characteristic 
length given by the width of the interface region. The scaling predictions are tested in 
the framework of the mean-field Ginzburg-Landau theory. Then we consider the Ising 
quantum chain in a linearly varying transverse field which corresponds to the extreme 
anisotropic limit of a classical two-dimensional Ising model. The quantum Hamiltonian 
can be diagonalized exactly in the scaling limit where the eigenvalue problem is the 
same as for the quantum harmonic oscillator. The energy density, the magnetization 
profile and the two-point correlation function are studied either analytically or by exact 
numerical calculations. Their scaling behaviour are in agreement with the predictions 
of the scaling theory. 



1. Introduction 

Inhomogeneities may have a strong influence on the properties of a system in the vicinity 
of a second-order phase transition. Actually this influence will depend on the relevance 
of the perturbation introduced by the inhomogeneity (see p] for a review). A relevant 
inhomogeneity may change the universality class of the system or even suppress the 
critical point as, for instance, in a finite-size system [2]. The critical behaviour may 
be altered by the presence of quenched disorder [3] or aperiodic modulation of the 
couplings [HE]. It will be modified locally (i.e., within a correlation length) at a flat 
free surface, at a corner [6-12] or at the tip of a parabolic-shape system [13], [T]. Line 
defects may have also some influence on the local critical behaviour [14-16,1]. One may 
mention the case of films where the presence of boundaries, by breaking translation 
invariance, leads to the formation of specific profiles [17-24]. 

Other types of inhomogeneities are linked to the application of external fields like 
magnetic, gravitational or thermal fields. These fields influence the behaviour of physical 
quantities like the magnetization, the particle or the energy density in the vicinity of the 
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homogeneous system critical point. One should first mention an early work on the xy- 
quantum chain in a linearly varying z-field [25] for which exact results were obtained in 
different scaling limits. Phase coexistence induced by antiparallel magnetic fields at the 
surfaces of a film was studied in [26-30]. The effect of gravity was considered in [28-30] 
where it was found that it restores two-phase coexistence up to the bulk critical point, 
above the wetting temperature. The effect of a temperature gradient on an interface 
at equilibrium was considered in [31] for a symmetric binary system below its critical 
point. Phase separation induced by temperature gradients was studied in [32-34]. The 
effect of temperature gradients on interfacial premelting was considered in [35] . 

The influence of inhomogeneities on the critical behaviour was also considered in a 
series of works on gradient percolation [36-40] where the perturbation was introduced 
as a tool allowing for high-precision estimates of the percolation threshold and the 
percolation exponents. 

In this work we begin with a presentation of the scaling theory for the density 
profiles in the presence of gradient field inhomogeneities. Specifically, we consider a 
system with a deviation from the critical coupling which varies linearly in one space 
direction. The coupling is at its critical value in the middle of the system where an 
interface region separates the ordered phase on the left from the disordered phase on 
the right. We test the validity of the scaling arguments, first at the mean-field level, 
within Ginsburg-Landau theory. Then we present a study of the Ising quantum chain 
in a linearly varying transverse field, hi — 1 + gl, which corresponds to the extreme 
anisotropic limit of the two-dimensional classical Ising model with a linear variation 
of the couplings. We work in the scaling limit where the size L of the system goes 
to infinity while the gradient g goes to zero with the product gL held fixed. The 
excitation spectrum of the inhomogeneous Ising quantum chain is obtained exactly in 
terms of the solution of an harmonic oscillator eigenvalue problem. The knowledge of the 
eigenvectors allows us to obtain the energy density profile, the magnetization profile and 
the behaviour of the spin-spin correlation function. Their scaling forms are in complete 
agreement with the results of the scaling theory. 

The paper is organized as follows: in section 2 we present a scaling analysis for 
the density profiles which is confirmed in section 3 in the framework of a mean-field 
approximation. Section 4 deals with the study the Ising quantum in a linearly varying 
transverse field. First the form of the quantum Hamiltonian is deduced from the 
classical Ising problem. Then the Hamiltonian is diagonalized exactly in the scaling 
limit. In the next subsections, analytical and numerical results for the energy density, 
the magnetization profile and the correlation function are confronted to the results of 
the scaling theory. We summarize our results in the last section. 
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2. Scaling arguments 

Let us consider a critical system perturbed by a constant gradient g along the ^-direction, 
such that the deviation from the critical coupling is given by 

K(z)-K c = A(z)=gz, g>0. (2.1) 

Thus around the origin there is an interface between the ordered phase on the left side 
where A < and the disordered phase on the right side where A > 0. Let £ denote the 
width of this interface. If the system is infinite £ is expected to depend on g only and 
to diverge when g vanishes. 

Under a change of the length scale by a factor b, the thermal perturbation A (z) 
with scaling dimension y t = 1/V, where v is the correlation length exponent, transforms 
as 



9'z = b^gz = </| , (2.2) 



so that 



9 = b 1+l/v g . (2.3) 
The interface width transforms as 

£' = £{g>) = - h =t (b 1+1 l»g) , (2.4) 

with b = g~ v '^ 1+v \ one finally obtains 

£ oc g-^'^ , (2.5) 

for the typical length introduced by the thermal gradient g. 

The same result can be obtained self-consistantly [36] by noticing that to the width 
£ is associated a typical deviation from the critical coupling A(£) = gi from which a 
characteristic length [A(£)]~ v can be constructed. Since the only length in the problem 
is the interface width, £, it satisfies 

t<x{gt)- v . (2.6) 

Solving for £, one immediately recovers (12.51) . 

Let us now study the influence of the thermal gradient on the scaling behaviour of 
the density (p with scaling dimension x v . This density can be the magnetization density 
m or the singular part of the energy density e. In the perturbed critical system it is a 
function ip(z,£) or, alternatively, according to (12. 5p . a function ip(z,g) transforming as 

<// = <pV,gr)=V+<p(z,g) (2.7) 
under a change of scale, so that 

^g) = b- x ^(^-,b 1+l/u g) ■ (2.8) 
With b = g~ u /( 1+u ) oc £ one obtains the scaling form 

<p(z,g) = g vx */( 1+v )$ \g v /( 1+v h] . (2.9) 



Gradient critical phenomena in the Ising quantum chain 



4 



One may notice that, according to (12. 5p . the prefactor in (12.91) exhibits the finite-size 
behaviour <p oc £~ x<p expected for a critical system with transverse size £. 

For a magnetic system the magnetization density, m(z,g), with scaling dimension 
x m — fl/V) is non- vanishing in the ordered region z < and one expects the same local 
critical behaviour as for the homogeneous system at the same value of the coupling, i.e., 

m{z,g) oc \A(z)f oc \gz\ UXm , z<0. (2.10) 

for not too large values of | A(z) | . Here and in the following, we suppose that m(z, g) > 0. 
Comparing fl 2 . 1 [) to (12.91) with (p — m, one obtains the form of the scaling function in 
this region: 

$ m (w) ~ \u\ UXm , u<0. (2.11) 

Just at the interface one expects $ m (0) = const. 7^ so that the local magnetization is 
non- vanishing and behaves as the prefactor in (12.91) . displaying the finite-size behaviour 
m(0,g) oc £~ Xm . Since m(z,g) > m(0,g) when z < 0, the scaling behaviour given 
in equations (12. lOj) and (12. lip is valid only when \gz\ UXm > £~ Xm . Using (I2.5p . this 
translates into z < —£. 

On the right of the interface, in the weak-coupling region, the magnetization density 
is expected to display the same exponential decay as the two-point correlation function 
r(z) = (m(O)m(z)). For a homogeneous system in its disordered phase, with a constant 
deviation A from the critical coupling, one has 

T(z) ~ exp ^— const. ~ exp ^— const. > z ^ £ > (2-12) 

where £ is the correlation length. Assuming the same behaviour for the inhomogeneous 
system, with the correlation length taking a value governed by the local deviation from 
the critical point, £ oc [A(,z)]~ I/ , one obtains [13] 

m(z, g) ~ exp (-const. g v z l+v ^j . (2.13) 

Equations (I2.9p and (12.131) give the form of the scaling function for the order parameter 
when z 3> 

m {z, g) = g^l^)^ ^/^z] , $ m (n) ~ exp (-const. u l+v ) . (2.14) 

The same type of arguments leads to the scaling behaviour of the singular part of 
the energy density 

e( z ,g) = ^/( 1 +-)cb c [g'/^z] , $ e (u) ~ exp (-const. \u\ 1+,/ ) , (2.15) 

for both sides of the interface. Here x e = d— 1/v is the scaling dimension of the energy 
density. 

The two-point correlation function, with scaling dimension 2x m , can also be written 
under the scaling form 

T(z,g) = g 2vx ^^ +v ^ v [g u/{1+u) z] , (2.16) 

with the exponential decay, $r(«) ~ exp (—const. u 1+u ), for the connected part, 
T c (z,g) = T(z,g) - {m(0))(m(z)). 
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3. Mean-field theory 

In order to check our scaling assumptions, let us now study the gradient perturbation 
problem in mean-fied theory. With a scalar order parameter m and up-down symmetry, 
the Ginzburg-Landau free energy functional of the critical system perturbed by the 
gradient term reads 



G[m] = G[0) 



v 



C,„ , 2 A{z) 2 B 4 ' 

— VmV H —m 2 H m 4 - ffm 

2 V ; 2 4 



tzy, (3.i) 



where B and C are positive constants. The first term is the energy contribution coming 
from inhomogeneities. The quadratic term is the thermal gradient perturbation with, 
as before, A(z) = gz and g > 0. The quartic term ensures the stability of the system in 
the ordered region for z < 0. The last term gives the interaction with the external field 
H. 

In the mean-field approximation the equilibrium value of the order parameter 
minimizes G[m]. Thus the variation of the free energy, 5G, vanishes to first order 
in 5m. This leads to the Ginzburg-Landau equation 

-C— r +gzm(z)+Bm 3 (z)=0 J (3.2) 
dz z 

where it was assumed that H = and that translation invariance is broken only in the 
z-direction. 

Introducing the dimensionless variable, ( = z/£, the Ginzburg-Landau equation 
can be rewritten as 

" + ^Cm(C) + Bm\Q = . (3.3) 

The coefficients of d 2 m/d( 2 and m(Q have to scale in the same way. Thus C/£ 2 cx gl 
and the interface width has the following scaling behaviour 

£~ (jj 1 3 cx g-»'^ , (3.4) 

since v — 1/2 in mean-field theory. 

In the ordered region, for z — £ or £ — 1 in (13. 3p . one can neglect the 
second derivative which is much smaller than the next term in the Ginzburg-Landau 
equation (13. 2p so that we obtain 

, . f-gz\ 1 / 2 , 
m z >9) ^ ( -g- J , (3.5) 

which confirms the scaling prediction (12.101) since here (3 = z/x m = 1/2. 

The order parameter is very small for z ^> and the cubic term in (13.21) can be 
neglected. With the change of variable u = z(gfC) 1 ^ 3 one obtains the Airy equation 

d 2 m 

—— — um(u) = . (3.6) 
du 2 
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The order parameter is then given by m(z,g) oc Ai(tt) with the leading asymptotic 
behaviour 

m(z, g) ~ exp const. u 3 ^ 2 ^j ~ exp(— const. g l ^ 2 z 3 ^ 2 ) , (3.7) 
when z 3> in agreement with equation (I2.14p with v = 1/2. 

4. Ising quantum chain in a linearly varying transverse field 

4-1. Quantum Hamiltonian 

Let us consider the two-dimensional classical nearest-neighbour Ising model on a square 
lattice with vertical couplings K\(l) varying in the horizontal direction {—L/2 < I < 
L/2) and constant horizontal couplings K 2 . The partition function can be written as 

Z oc TrT M , (4.1) 

where M is the number of horizontal rows and T is the row-to-row transfer matrix given 
by HUES] 



T = exp 



exp 



(4.2) 



E^X+i 
. i 

The cr 2 ^ are Pauli spin operators, A*(/) = — 1/2 lntanh[Ai(Z)] is the dual of the vertical 
coupling and it is assumed to vary as 

K* 1 {l)=K* 1 {l + gl), g>0. (4.3) 

In the extreme anisotropic limit [43-45] K\ — > oo (so that K* — > 0) and A 2 — > 
while keeping the ratio fa = KI/K2 constant, the transfert matrix T can be rewritten 

as 

T = 1 - 2K 2 H . (4.4) 

7^ is the Hamiltonian of the Ising quantum chain in a transverse field jl6] and takes the 
following form 

L/2-1 1 L/2 

n = -- E "fof+i-, E ^ = mi + ^), (4.5) 

Z J=-X/2 / l=-L/2 

where g = 9/L with 9 > 0. 

In the homogeneous case, g = 0, the system is self-dual and critical at fa = 1 in the 
thermodynamic limit L — > 00. For fa < 1 the system is ordered with (o~ x ) 7^ whereas 
(a x ) = when fa > 1. 

The perturbed system is assumed to be critical when 9 = so that fa = 1 in (14.51) 
and hi = 1+91/ L = 1+gl. In order to keep fa/ > 0, one has to take 9 < 2. The transverse 
field is smaller (greater) than its critical value on the left (right) of the origin. Thus, in 
the thermodynamic limit, the quantum chain is ordered for I < and the value of the 
order parameter decreases to zero when I — > +00. 
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4-2. Diagonalization 

After a Jordan- Wigner transformation [47] , the Hamiltonian ( 14.51) becomes a quadratic 
form in fermion creation and annihilation operators which is diagonalized through a 
canonical transformation [46j HH] leading to 

n = jZe q {n\n q -\). (4.6) 

<?=0 A 

The 77T (7^ ) are diagonal fermion creation (anihilation) operators and the e q are the 
energies of the fermionic excitations. They satisfy the following set of equations 



-qYq 1 



A f g = £q*/>, 



where 



/ -h. 







L/2 





-h. 






qyq 1 


L/2+1 

1 





(4.7) 













-h L /2-i 



-h 



(4- 



L/2 / 



and is the transposed matrix. According to (14. 7p . the normalized eigenvectors, <p q 
and ipq, are solutions of the following eigenvalue equations 



aaV„ 



-q^q > 



(4.9) 



In the bulk of the system these eigenvalue equations can be written as 



- 1) + 



hf <j> q {l) + hi<t> q {l 



h t ij) q (l-l)+ e 2 -l-h 2 i) q (l)+h l+1 ip q (l + l) 



0. 

: 0. 



(4.10) 



Introducing the scaling variable u = g 1//2 /, the difference equations in (I4.10p can be 
expanded in the scaling limit where L — > 00 and g — > in such a way that the product, 
gL = 0, is held fixed. Up to terms of the first order in g the expansions lead to the 
following harmonic oscillator eigenvalue equations 



d 2 



dv/ 



+ 



JJ 



1/2 



- 1 



u 



(f)(u) = . 



dSp_ 

du 2 



+ 



JJ 



1/2 



+ 1-U 2 



ip(u)=0. (4.11) 



With the boundary conditions 0(±oo) = ^>(±oo) = 0, one obtains 

^Ju) = C n e- u2 / 2 H n (u) , 



4> n +i(u) = ip n (u) , e n = ^/2ng, n = 0,1,2,..., (4.12) 

where H n (u) is the Hermite polynomial of order n and C n = (g/-K) 1 ^ 4: (2 n n\)^ 1 ^ 2 is a 
normalization factor. 

One may notice that in the bulk, according to (14.81) . the operators A and A^ satisfy 
the commutation relation: 



[A, A 1 "] = gC , C ltm = S, 



l— l,m 



+ 5, 



Z+l,m 



(4.13) 
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Introducing the normalized operators, a = A/y/2g and = A^/y/2~g, in (14.131) and 
applying the commutator to a test function, in the scaling limit one obtains the canonical 
bosonic commutation relation 

[a,at] = l + 0(<7). (4.14) 

Thus, once normalized in the same way, the equations in (14. 7p correspond to the lowering 
and raising operations on the eigenstates ip n and <p n of the harmonic oscillator. 

The eigenstate (f>o(u), associated with the excitation e which is vanishing in the 
scaling limit, cannot be obtained from the harmonic oscillator equation since it is 
incompatible with the boundary condition at —00. It corresponds to a mode which is 
localized in the vicinity of the left boundary and which is related to the presence of a non- 
vanishing magnetization m s at the left boundary of the system [501 [51]. In order to obtain 
the form of this localized mode and the actual value of the corresponding excitation, 
one has to come back to the original finite-size system with open boundary conditions. 
It can be shown that, when it vanihishes faster than 1/L, the lowest excitation is given 
by [52-54] 

e ~ m s m* , (4.15) 

where m* is the magnetization of the dual chain at the left boundary. m s and m* are 
equal to the component (f> (—L/2) of the normalized eigenvector and can be deduced 
from the second equation in (14.71) with e — 0. These boundary magnetizations take the 
following forms [50] 



L/2 k 

1+ e n 

k=-L/2 l=-L/2 



-1/2 



L/2 k 

1+ e n k 2 

k=-L/2l=-L/2 



-1/2 



(4.16) 



With hi = 1 + gl it is straightforward to show that m s is 0(1) whereas the dual 
magnetization is exponentially small: 

m s ~ I — exp — . (4-17) 



, 7T , 

As a consequence, the smallest excitation vanishes exponentially with the size of the 
system as 

e ^e' 9L2/s . (4.18) 

Using the second equation in (I4.13P with e ~ 0, one obtains 

M-L/2 + I) ~ M-L/2)e- l ' l ° , lo 1 = | ln(l - gL/2)\ . (4.19) 

The localisation length Iq diverges as (6/2)~ l when the local deviation from the critical 
transverse field, 9/2 = 1 — h-L/2, vanishes. It behaves as the local correlation length 
since v = 1 for the Ising quantum chain. In the scaling limit g — > and L — > 00 while 
9 = gL remains constant; thus the localization length l remains also constant and the 
localized mode, (f>o{u), has a vanishing amplitude in the bulk. 

This completes the diagonalization of the model in the continuum limit. 
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Figure 1. Rescaled energy density profiles for different chain sizes L and 9 values. 
The corresponding values of the gradient are g = 1 x 1 0~ 3 , 5 x 1 0~ 4 , 2 x f 0~ 4 , f x f 0~ 4 
from top to bottom in the legend. The solid line is the analytic result in (|4.24|) . 



4-3. Energy density profile 

Let us consider the connected autocorrelation function in imaginary time 

G c (l,r) = (af(r)af)-(af) 2 , (4.20) 

where erf (r) = e rn afe~ Tn . It can be rewritten as the eigenstate expansion 

Ge(/,r)=EK^f|0)|V^-^, (4.21) 

where |0) is the ground state of 7i with eigenvalue E and \i) an excited state with 
eigenvalue E^. Since the operator erf is a two-fermion operator, the only non- vanishing 
matrix elements in the expansion are the two-fermion states |z) = T^ry^S). When r — > oo, 
the amplitude of the dominant term defines the off-diagonal energy density 

e(0 = Ke|af|0>|, (4.22) 

where |e) = ?7i?7o|0) is the first even excited state. The expansion of erf in terms of 
diagonal fermions leads to |56j 

e(/) = |^i(/)0 o (O-^o(O0i(OI- (4-23) 
Inserting (I4.12p into (I4.23p . one finally obtains the Gaussian profile 

e (/) = |^o(0| 2 = /fexp(-^ 2 ), (4.24) 

which is in agreement with the form (12.151) deduced from scaling considerations since 
v = x c = 1 in the Ising model. 

In order to reduce the boundary effects due to the localized mode, one has to 
consider sizes L ^> l ~ O^ 1 for small 9. Typically, for 9 = 0.1, Iq ~ 20 and taking 
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L > 10 x l is sufficient to suppress the localized mode effect. The behaviour of the off- 
diagonal energy density, deduced from (I4.23P through numerically exact diagonalization 
for chains of size up to L = 1000, is shown in figure [U An excellent agreement with the 
analytic result in (14.241) is obtained. 
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Figure 2. Rescaled magnetization profiles for different chain sizes L and 9 values. 
The values of the gradient are the same as in figure Q] 



4-4- Magnetization profile 

Due to the symmetry of the Hamiltonian, the ground state expectation value of 
of identically vanishes. In order to obtain the magnetization profile one has to break 
this symmetry. Since the ground state |0) is asymptotically degenerate with the lowest 
one-fermion state \a) = 770" 1 0) such that E a — E = e — 0, one can construct the state 

l+) = ^(|0> + k» (4.25) 

and take the expectation value of erf in that state which leads tdj] 

m(l) = <+kf|+> = (a\af\0) . (4.26) 

One may also notice that the imaginary time autocorrelation function has the following 
asymptotic behaviour 

fim G m {l,r) =m 2 (l) = (cr\af\0) 2 , (4.27) 

which is another way to justify the use of the off-diagonal matrix element for the 
magnetization profile. 

I More generally, with the linear combination la) = l°)+ a l cr ) , where a = pe l8 , one obtains (a\af \a) = 
2 i1°i^ e ( a \ a i |Q) which is maximum for a = 1 when (a\af\ijj a ) = (a\af\0). 
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Figure 3. Behaviour of the magnetization profile for z < 0. The dashed line gives 
m 8 = —2gz, the behaviour expected from equation (|2.10p . valid in the scaling limit 
where g — > 0. The solid line corresponds to to 8 = — 2gz — g 2 z 2 (see text). The values 
of the gradient are the same as in figure [1] 




gz 

Figure 4. Semi- logarithmic plot of the rescaled magnetization profile in the disordered 
region, z > 0. A linear behaviour is expected from equation (|2. 14|) . The values of the 
gradient are the same as in figure [TJ 



Rewriting af in terms of diagonal fermions and using Wick's theorem, the local 
magnetization can be expressed as a determinant [55] 



H 



L/2 
L/2+1 



Hi 



L/2 -L/2 
L/2+1 -L/2 



Gl-l,-L/2 
Gl-L/2 



L/2, -L/2+1 



L/2+1, -L/2+1 



Gl-l- 



L/2+1 



Gl -L/2+l 



L/2,1-1 
L/2+1, Z-l 



Gl-ij-i 

Gij-i 



(4.28) 
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where 

H, = MJ) , Gj, k = -E Mj)Mk) ■ (4.29) 

n=0 

Figure [2] shows the rescaled magnetization profile obtained for chain sizes up to L = 1000 
and different values of 9. The numerical results are in excellent agreement with the 
scaling behaviour of equation (12.141) with v — 1 and x m = 1/8 for the Ising quantum 
chain in a transverse field. 

In order to test the scaling assumption f)2.10p . m 1 ^ ram ^ ) = m 8 is plotted as a 
function of gz in figure [3j The spontaneous magnetization of the homogeneous Ising 
chain varies as m = (1 — h 2 ) 1 ^ 8 as a function of the transverse field [46J for h < 1. 
Replacing h by its local value h(z) = 1 + gz in the inhomogeneous system, one obtains 
[m(z, g)] 8 = —2gz — g 2 z 2 for z < in very good agreement with the numerical data. 
One may notice that the second term in this expression does not fit with the scaling 
form which follows from ( 12. 9p but this term becomes negligible in the scaling limit where 
g — > as may be verified on figure [3j 

Figure H] shows the behaviour of the rescaled magnetization profile in the 
disordered region on a linear-logarithmic scale in order to check the exponential decay, 
exp (— const. gz 2 ), expected from equation ( 12.141) . The exact numerical results are once 
again in quite good agreement with the expected scaling behaviour. 

4-5. Correlation functions 




-2-10 1 2 3 

1/2 

8 Z 

Figure 5. Scaling behaviour of the spin-spin correlation function T(z, g) for g = 
5 x 10~ 4 (diamond), 1 x 10~ 4 (square). The dashed line corresponds to the product 
(m(Q))(m(z)) and the inset gives the behaviour of the connected part. 
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Figure 6. Semi-logarithmic plot of the rescaled connected part of the spin-spin 
correlation function as a function of gz 2 . The linear behaviour is in agreement 
with [[432]) , 



The correlation between a spin at / and the central spin is measured by the 
correlation function 



r(j) = (okxio) 

given by the determinant [46] 

r(0 



(4.30) 



^2,0 G 



2,1 



GYo Gi 



G\,i-\ 
G^i-i 

Gu-i 



(4.31) 



which involves the contractions already defined in (14.29ft . 

The numerical results, shown in figure [51 confirm the scaling form given in (I2.16p . 
The behaviour of the connected part is shown in the inset. A closer analysis of the decay 
of T c (z,g) points to the following Gaussian behaviour 



-1/4, 



exp I —const. gz' 



(4.32) 



in agreement with the scaling prediction ( 12. 16ft . The constant takes a different value on 
the two sides of the system (see figure [6]). 



5. Summary and conclusion 

We have presented a scaling theory for the behaviour of the magnetization profile, the 
energy density profile and the two-point correlation function in a critical system in 
the presence of a linearly varying deviation from the critical coupling, A(z) = gz. The 
gradient g introduces a new length scale in the problem, I ~ g^ u /( 1+u )^ which depends on 
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the correlation length exponent v in the case of a thermal perturbation considered here. 
The form of the scaling functions have been obtained by assuming that the different 
physical quantities have locally the same functional form as in the homogeneous system 
with the deviation from the critical point, A, replaced by its local value, A(z), in the 
inhomogeneous system. 

The results of the scaling theory have been confirmed, first in mean-field theory 
and then in a study of the Ising quantum chain in a linearly varying transverse field. In 
this later case the excitation spectrum of the quantum Hamiltonian has been obtained 
exactly in the scaling limit where the size of the system L — > oo, the gradient g — > 
while the product Lg is held fixed. In this continuum limit one recovers the eigenvalue 
equation of the harmonic oscillator problem. The knowledge of the eigenvectors allows 
us to calulate exactly the energy density profile and to obtain numerically exact results 
for the magnetization density and the two-point correlation function. The very good 
agreement with the scaling results strongly supports the locality assumption used to 
deduce the form of the scaling functions. 
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